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ABSTRACT 

We present radiometric effective radii and visual geometric albedos for six asteroids 
in comet-like orbits. Our sample has three of the four known retrograde asteroids (1999 
LE31, 2000 DCs, and 2000 HE46) and three objects ((18916) 2000 OG44, 2000 PG3, 
and 2000 SBi) on prograde but highly elliptical orbits. These measurements more than 
double the number of known albedos for asteroids with a Tisserand invariant in the 
cometary regime. We find that all six of our objects, and nine of the ten now known, 
have albedos that are as low as those of active cometary nuclei, which is consistent 
with their supposed evolutionary connection to that group. This albedo distribution is 
distinct from that of the whole near-Earth and unusual asteroid population, and the 
strong correlation between Tisserand invariant and albedo suggests there is a significant 
cometary contribution to this asteroid population. 

Subject headings: asteroids — comets: general 



Introduction 



An old comet that loses all of its available volatiles or is covered by a mantle that prevents 
sublimation of subsurface ice will observationally appear as a near-Earth or unusual ^ asteroid 
(NEA or UA). The dynamical lifetime of short-period comets are about 10 to 100 times longer 
than the devolatilization timescale (Levison & Duncan 1994), so one expects to see such extinct 
comets if they do not disintegrate or collide with a planet. However there is currently no way to 
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in transjovian space. 
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determine if a given asteroid is such a comet (short of witnessing last gasp activity) , so we do not 
yet know the fraction of extinct comets in the asteroid population. Knowing this fraction would 
give clues to the physical evolution of comets and the hazard to Earth from asteroid and comet 
collisions. 

Dynamical models of the cometary component are rendered inconclusive by the unknown extent 
and duration of nongravitational forces on active nuclei (Harris &: Bailey 1998). Some models 
(Bottke et al. 2000) can even fit the NEA population without a cometary source. Nevertheless, 
the existence of the now-inactive comet 107P /Wilson-Harrington, which has been observed to have 
a coma just once in the 51 years since discovery (Bowell 1992; Fernandez et al. 1997), strongly 
suggests that there is a non-zero cometary contribution, and we have taken an observational tack 
to address the problem. A statistical indicator for the cometary origin of an asteroid comes from 
the albedo and orbit shape. All known albedos of cometary nuclei are low and almost all inner 
Solar System cometary orbits have common characteristics. Thus an asteroid in a comet-like orbit 
with a comet-like albedo is a good candidate extinct comet. 

Nearly a dozen albedos of active cometary nuclei have been measured and all are dark, with 
geometric albedos ranging from 0.02 to 0.12 (Jewitt 1992; Fernandez 1999). Thus we expect to see 
similar albedos among the extinct comet fraction of NEAs and UAs. Albedo alone is not a unique 

determinant of origin however since there are many non-comctary, dark asteroids in the outer Main 
Belt (Gradie, Chapman, k, Tedesco 1989) from which some NEAs and UAs could have come. 

The Tisscrand invariant Tj (Tisserand 1896), a constant of motion in the restricted three- 
body problem, can be used to separate objects by dynamical class. The threshold Tj = 3 separates 
objects coupled to (< 3) or decoupled from (> 3) Jupiter. Generally asteroids have Tj > 3, Jupiter- 
family (JF) comets have 2 < Tj < 3, and Halley-family (HF) and long-period (LP) comets have 
Tj <2 (Levison 1996), though the scheme is not fail-safe. There are currently (mid- April 2001) 11 
asteroids with Tj < 2, A of which are in retrograde orbits, and 131 NEAs and UAs with 2 < Tj < 3. 
The total number of NEAs and UAs is currently 1,400, so 10% have Tj < 3. Among NEAs alone, 
82 out of 1,327, or 6%, have Tj < 3. In this Letter we describe the radiometric determination 
of the albedos and effective radii of six Tj < 3 asteroids. The Tisserand values, the contributing 
orbital elements, and the geometry of our observations are given in Table 1. The new data bring 
to ten the total number of Tj < 3 asteroids with known albedos. 

2. Observations and Reduction 

The observations span two wavelength regimes, mid-infrared (MIR) and visible (simultaneously 
for 4 of the 6 objects). The MIR images were obtained with the Keck I telescope using the LWS 
array (Jones & Puetter 1993) in June 2000 and with the Keck II telescope using the MIRLIN array 
(Ressler et al. 1994) in November 2000. The visible data were obtained with the Univ. of Hawaii 
2.2-m telescope using a Tek2048 CCD in July and November 2000. Table 2 gives the averages of 
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the measured flux densities. All objects were point sources. 

The MIR data were obtained using chopping and nodding, with throws of 4 arcsec. Non- 
sidereal guiding was used for each target. Flat fields were obtained by comparing staring images 
taken at both high and low airmass. For photometric calibration of LWS data we compared count 
rates to the following known (12.5 and 17.9 /um) flux densities of standard stars: a Lyr, 26.4 Jy 
and 12.9 Jy; a Lib, 120.7 Jy and 58.9 Jy; a CrB, 3.64 Jy and 1.97 Jy; 7 Aql, 54.3 Jy and 27.5 Jy. 
For our MIRLIN data we used the following known (11.7, 12.5, and 20.8 iim) flux densities and 
stars: 7 Aql, 61.7 Jy, 54.3 Jy, and 20.6 Jy; /? Peg, 313 Jy, 277 Jy, and 107 Jy; a Ari, 62.9 Jy, 55.5 
Jy, and 21.1 Jy; a CMi, 59.7 Jy, 52.2 Jy, and 18.8 Jy. The values are derived from the standard 
system of Tokunaga (1984). Comparing raw photometry over a range of airmasses let us find the 
extinction corrections: 0.12 and 0.35 mag/airmass at 12.5 and 17.9 ;um respectively for the June 
data; 0.08, 0.10, and 0.40 mag/airmass at 11.7, 12.5, and 20.8 //m respectively for the November 
data. To maximize the signal-to-noise ratio in the photometry we used aperture corrections derived 
from nearby standard star radial profiles. 

The visible images were obtained while guiding on a nearby star with sidereal tracking rates 
in July but non-sidereal rates in November. A flat field was obtained by combining images of the 
blank twilight sky. Flux calibration and airmass corrections were calculated by measurements of 
the Landolt (1992) standard stars SA 107-457, -456 and the PG 1323-086 group in July 2000, and 
stars SA 98-966, -1002, -L3, and -L4 in November 2000. 



The basic radiometric method to obtain an effective radius R and geometric albedo p is to 
solve two equations with these two unknowns, first done by Allen (1970) and described in detail by 
Lebofsky & Spencer (1989): 



where F is the measured flux density of the object at wavelength A in the visible ("vis") or 
mid-infrared ("mir"); Fq is the flux density of the Sun at Earth as a function of wavelength; r 
and A are the object's heliocentric and geocentric distances, respectively; $ is the phase function 
in each regime; Bx is the Planck function; e is the infrared emissivity; 77 is a factor to account 
for infrared beaming; and T is the temperature, which is a function of p, surface planetographic 
position and the phase integral q which links the geometric and Bond albedos. For lack of 



3. Analysis 
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detailed shape and rotational information - as is the case for our six objects - the modeled body 
is assumed to be a sphere. 

The temperature is calculated using a model of the thermal behavior. Unfortunately, the 
thermal inertias are largely unknown so we use two widely-employed simple models that cover the 
extremes of thermal behavior: one for slow-rotators and one for fast. The former (a.k.a. "standard 
thermal model" , STM) applies if the rotation is so slow (or the thermal inertia so low) that every 
point on the surface is in instantaneous equilibrium with the impinging solar radiation. The latter 
(a.k.a. "isothermal latitude model" , ILM) applies if the rotation is so fast (or the thermal inertia so 
high) that a surface element does not appreciably cool as it spins away from local noon and out of 
sunlight. The extreme case occurs when the object's rotation axis is normal to the Sun-object-Earth 
plane. 

There are other parameters to the models: e, 77, ^mirj ^visj ^ind q. Emissivity is close to unity 
and we will assume a constant value of 0.9 here. The beaming parameter is unity by definition in the 
ILM. In the STM, the value is known for only a few asteroids and can range from approximately 0.7 
to 1.2 (Harris 1998; Harris &; Davies 1999). For our applications of the STM we assume a possible 
range of 0.75 < r/ < 1.25, uniformly distributed. 

For "I> we assume that the magnitude scales with the phase angle a: — 2.51og$ = /3a. In the 
MIR, the effect is only loosely constrained but based on earlier work with the STM (Lebofsky et 
al. 1986) we assume a range of 0.005 mag/deg < l3mir < 0.017 mag/deg, uniformly distributed. 
For the ILM f3mir = by definition. In the visible, studied in much more detail (e.g. by Lumme 
& Bowell (1981)), dark asteroids follow /3,,,f,5 ~ 0.035 mag/deg, and we assume here 0.025 mag/deg 
^ Pvis ^ 0.045 mag/deg. This coefficient also determines q, but since that has a minor effect on the 
modeling we will leave that parameter a constant 0.5, the integral's value for Pyis = 0.034 (Allen 
1973). 

Our radius and albedo results from this modeling are shown in Table 3. Since we have as 
many or more parameters in the model than data points, useful x^-distribution calculations are 
impossible. Instead, we found the range of valid R and p by sampling parameter space and declaring 

a good fit when each model flux density passed within 1.5-cj of its data point. The "error bars" in 
Table 3 actually show the full range (not 1-a limits) of values that yield acceptable fits. For 1999 
LE31, we used the average of the modeling results from the two nights of MIR photometry. 

The STM and ILM predict very different color temperatures, providing a way to differentiate 
them. For 2000 DGs, 2000 OG44, and 2000 SBi, only the STM provides acceptable fits. For 2000 
PG3, the ILM generally gives poorer fits and is acceptable only for a small range of radii and 
albedos. For 1999 LE31 and 2000 HE46, we only have one MIR wavelength and thus cannot find 
the color temperature. We note that the ILM albedos are very small, smaller than for any other 
Solar System object. This, combined with the generally poorer fits, implies that the objects are 
better interpreted as being slow rotators, although the heliocentric distance of 1999 LE31 may be 
high enough to make it a borderline case between the two extremes (Spencer, Lebofsky, & Sykes 
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1989). 

Some further caveats are notable. First, the values in Table 3 are valid in the context of the 
thermal models; the models represent extrema of thermal behavior and the ranges in Table 3 do 
not describe the systematic errors from the models themselves. However these errors are likely 
to be comparable to the formal errors so the values are physically meaningful. Second, rotational 
variation could corrupt the albedo calculations for the two asteroids with non-simultaneous MIR 
and visible photometry, which effectively increases the albedo errors. Third, our observations of 
2000 PG3 took place at a = 59°, which is beyond the range of any measured thermal phase function, 
so we do not know if the standard formula for ^mir mimics reality at such angles. 

The radii of our six objects are not unusual in comparison to cometary nuclei (Meech, Hainaut, 
&L Marsden 2000) and other, dynamically-similar asteroids (Veeder et al. 1989), and we now turn 

our attention to the albedos. Figure 1 displays the albedos and Tisserand invariants of our six 
objects as well as all other, 42 known albedos of NEAs and UAs. Also plotted are albedos of active 
cometary nuclei. 

Whereas previously there were only four Tj < 3 asteroids on the plot, there are now ten, 
and all six new ones fall squarely within the cometary regime of albedos. Moreover, with these 
new additions we now see evidence of a trend with Tj. Ninety percent of the Tj < 3 asteroids 
have comet-like albedos (median albedo is 0.038), whereas only two of the 38 Tj > 3 objects do 
(median albedo is 0.215). An important bias in the distribution here is that it is easier to discover 
shinier asteroids and thus they are overrepresented in the sample, but nevertheless the difference 
in distributions is quite stark. 

The break in the distributions occurs near the line Tj = 3, the cutoff for Jupiter coupling, which 
suggests a dynamical relationship with this effect. The low albedos of nearly all Tj < 3 objects 
are consistent with a significant fraction of extinct comets among the NEA and UA population. 
Specifically, if 10% of NEAs and UAs have Tj < 3, and 90% of those have comet-like albedos, the 
extinct comet candidate fraction would be approximately 9%. This is a lower limit to the actual 
fraction of candidates because of the albedo bias mentioned above. 



4. Summary 

Using MIR and visible photometry, and employing the widely-used standard thermal model 
(STM) for slow rotators, we have derived new effective radii and geometric albedos for six asteroids 
in comet-like orbits; all six have Tj < 3. We find the following: 

• All six objects are dark, as dark as the albedo spread of cometary nuclei. The radii are also 
similar to those of active nuclei. This is consistent with a cometary origin, as if the asteroids were 
formerly active comets that lost all near-surface volatiles. 

• For four objects we have photometry at 2 or 3 MIR wavelengths, and the STM yields an 
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excellent description of the color temperature, better than the fast-rotator model. 

• Plotting all 48 known NEA and UA albedos - including our 6 new ones - vs. Tj shows a 
markedly sharp break, virtually a step function, at the line Tj = 3. 

• Eleven of the 48 objects have comet-like albedos: fully 90% (9 of 10) of the Tj < 3 objects, 
but only 5% (2 of 38) of the Tj > 3 objects. The median albedos p and their r.m.s. scatters are: 



• This disparity in median albedo suggests the fraction of extinct comets among NEAs and 
UAs is significant (at least 9% are candidates) and that enough cometary nuclei have sufficiently 
long physical lifetimes to survive devolatilization without disintegrating. 

It is clear that further surveys of asteroid albedos are necessary. Specifically, only 4 of the 10 
Tj < 3 objects are NEAs; more members of that group need to be sampled. Furthermore, with 
the recent explosion in asteroid discoveries and the ready availability of sensitive MIR detectors, a 
less biased sampling of albedos should be undertaken to obtain a confident estimate of the albedo 
distributions. 

We are indebted to Michael Ressler for allowing MIRLIN's use on Keck and to Varoujan 
Gorjian for instrument support. The operation of MIRLIN is supported by an award from NASA's 
Office of Space Science. We thank operators Joel Aycock, Meg Whittle, Wayne Wack, and John 
Dvorak for their assistance. We acknowledge the JPL SSD group for their very useful "Horizons" 
ephemeris program. We appreciate the help of Olivier Hainaut in improving this manuscript. This 
work was supported by grants to DCJ from NSF. 
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Table 1: Orbits and Observing Geometry 



Object 


a 


e 


i 


Tj 


UT Date 


r 


A 


a 




(AU) 




n 




(A.D.2000) 


(AU) 


(AU) 


n 


1999 LE31 


8.16 


0.472 


152 


-1.31 


Jun 22 


5.238 


5.118 


11.2 


n 










Jun 23 


5.240 


5.137 


11.2 


n 










Jul 2 


5.267 


5.311 


11.0 


2000 HE46 


24.3 


0.903 


158 


-1.51 


Jun 23 


2.526 


2.470 


23.5 






55 


55 


55 


Jul 2 


2.562 


2.689 


22.2 


2000 DGs 


10.8 


0.793 


129 


-0.62 


Nov 8 


2.319 


1.728 


22.9 


2000 OG44 


3.88 


0.581 


7.33 


2.74 


Nov 8 


1.665 


0.920 


30.6 


2000 PG3 


2.83 


0.859 


20.5 


2.55 


Nov 8 


1.065 


0.929 


59.1 


2000 SBi 


3.34 


0.541 


22.2 


2.81 


Nov 8 


1.554 


0.673 


25.6 



Note. — a = semimajor axis, e = eccentricity, i = inclination, r = heliocentric distance at time of observation, A 
= geocentric distance at time of observation, and a = phase angle at time of observation. 
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Table 2: Asteroid Pliotc^moiry 



Object 


UT Date 


UT Time 


Wavelength 


Flux Density 


Meas. 




(A.D.2000) 






(mjy, Jy, or mag) 


b 


1999 LE31 


Jun 22 


07:11-07:27 


12.5 


6.35 ± 0.78 mJy 


2 




Jun 23 


06:32-06:49 


12.5 


5.65 ± 0.71 mJy 


2 




Jun 23 


06:17 


17.9 


< 57 mjy'= 


1 




Jul 2 


07:32 


0.65 


20.44 ± 0.05 mag 


1 


2000 HE4fi 


Jun 23 


08:01-08:20 


12.5 


29.4 ± 1.9 mJy 


2 




Jul 2 


06:16 


0.65 


20.11 ± 0.02 mag 


1 


2000 DCs 


Nov 8 


13:06-13:12 


11.7 


0.382 ± 0.013 Jy 


4 




Nov 8 


12:35-13:38 


12.5 


0.411 ± 0.037 Jy 


13 




Nov 8 


12:43-13:28 


20.8 


0.71 ±0.15 Jy 


3 




Nov 8 


12:12-12:47 


0.65 


16.826 ± 0.016 mag 


8 


2000 OG44 


Nov 8 


06:07-06:12 


12.5 


0.766 ± 0.032 Jy 


3 




Nov 8 


06:14-06:20 


20.8 


0.739 ± 0.079 Jy 


4 




Nov 8 


06:15-06:17 


0.65 


16.39 ± 0.01 mag 


2 


2000 PG3 


Nov 8 


04:39-05:13 


12.5 


0.650 ± 0.094 Jy 


6 




Nov 8 


04:43-05:15 


20.8 


0.60 ± 0.12 Jy 


2 




Nov 8 


05:11-05:16 


0.65 


17.857 ± 0.013 mag 


2 


2000 SBi 


Nov 8 


11:07-11:09 


11.7 


1.112 ±0.035 Jy 


2 




Nov 8 


11:02-11:48 


12.5 


1.139 ±0.061 Jy 


8 




Nov 8 


11:11-11:40 


20.8 


1.48 ± 0.35 Jy 


3 




Nov 8 


11:15-11:18 


0.65 


16.27 ± 0.01 mag 


2 



""0.65" refers to R band. 

"Meas." signifies the number of measurements that were used to calculate the value in the "Flux Density" column. 
°This is a 3-a upper limit. 



- 11 - 



Table 3. Effective Radii and Geometric Albedos' 



Object 


Slow Rotator Model 


Fast Rotator Model 




Eff. Radius (km) 


Geom. Albedo 


Eff. Radius (km) 


Geom. Albedo 


1999 LE31 

2000 HE46 
2000 DGs 
2000 OG44 
2000 PG3 
2000 SBi 


9.05lt°t 

q cc+l.lO 

0.00_Q yg 

07+0.50 

<J-o'_0.40 
nS+1-42 
^■'-'°-0.95 
q r7+0.92 
-0.62 


Q QO. +0.030 

0.0271°:°?^ 

0.0381H}? 
021+° °^^ 

U.U±y_QQ^Q 


00 07+2.93 

D.OO_Q 28 

b 

b 

q 4q+0.21c 
^•^^-0.19 
b 


0041+°-°°^^ 
0067+° °°22 

U.UUDr_Q QQ21 
b 

b 

n n-i c:+0.007c 

U.U±0_Q 

b 



^With full ranges of acceptable values (not l-cr limits). 
''The fast model gives unacceptable fits for these objects. 

'^The fast model gives generally poorer fits and is acceptable for only a small spread 
of radii and albedos. 
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Fig. 1. — Plot of Tisserand invariant vs. all known geometric albedos for Tj < 3 NEAs and UAs 
(green circles), Tj > 3 NEAs and UAs (blue hourglasses), and comets (red squares). The six 
objects presented in this Letter arc marked as filled circles with estimated 1-cr errors. A heavy 
vertical line marks the dynamical boundary Tj = 3 (see §1). Of the 48 total asteroids plotted, 
11 have comet-like albedos, or 23%, but fully 90% (9 of 10) of those with Tj < 3 have comet-like 
albedos. This is consistent with a cometary origin for those asteroids and a significant cometary 
contribution to that dynamical group of asteroids. The plotted data were obtained from this work, 
Morrison, Gradie, & Rickc (1976), Cruikshank & Jones (1977), Lebofsky et al. (1978), Lebofsky et 
al. (1979), Lebofsky et al. (1981), Green, Meadows, & Davies (1985), Tcdesco & Gradie (1987), Bell, 
Hawke, & Brown (1988), Veeder et al. (1989), Tedesco (1992), Hudson & Ostro (1995), Mottola 
et al. (1997), Pravec et al. (1997), Harris (1998), Thomas et al. (2000), Harris & Davies (1999), 
Fernandez (1999), Delbo et al. (2000), and Fernandez et al. (2001). 



